When is Quantum Decoherence Dynamics Classical? 
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A direct classical analog of quantum decoherence is introduced. Similarities and differences be- 
tween decoherence dynamics examined quantum mechanically and classically are exposed via a 
second-order perturbative treatment and via a strong decoherence theory, showing a strong depen- 
dence on the nature of the system-environment coupling. For example, for the traditionally assumed 
linear coupling, the classical and quantum results are shown to be in exact agreement. 
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■ Decoherence is the loss of quantum coherence due to system-bath coupling. There has been considerable theoretical 
' and experimental work demonstrating that quantum-classical correspondence (QCC) can be induced by decoherence 
, 10, H, p. By contrast, little work has been done on examining the correspondence between classical and quantum 

■ descriptions of the time evolution of decoherence itself, i.e. decoherence dynamics. In this Letter we show that (a) one 
' can introduce a direct classical analog of quantum decoherence, and (b) by examining the dynamics of decoherence 

classically one gains new insights into both the dynamics of decoherence described quantum mechanically and into the 
conditions for QCC of the dynamics of decoherence. For example, we show that the extent of QCC depends strongly 
on the nature of the system-bath coupling and far less upon Ti than expected, that results assumed to be quantum 
mechanical can be obtained classically and that nonlinear system-bath coupling can cause nonclassical decoherence 
^ , dynamics even for macroscopic systems. 

The formal Liouville-based theory of QCC in an isolated system |Q makes clear that there is a strict analogy 
between quantum and classical dynamics in phase space. As in the quantum case, the classical Liouville dynamics of 
' a closed system is unitary and we expect the reduced classical Liouville dynamics of a system coupled to a bath to 
] be nonunitary. We therefore suspect that, due to bath's coarse-graining effects, the reduced dynamics of the system 
fS| propagated classically will show decoherence dynamics that is, qualitatively, parallel to that seen in quantum dynamics 
' insofar as the loss of phase information, entropy production, etc. 

, Here we quantitatively compare the dynamics of decoherence that is induced quantum mechanically to that induced 
^ ^ classically. This is done by analytically examining the dynamics of an initial quantum state, in a system coupled to a 
^ I ^ , bath, that is propagated either quantum mechanically or classically. The observed (rather remarkable) similarities and 
(~| ■ differences between the classical and quantum decoherence dynamics should be of considerable interest to a variety of 
modern fields such as quantum information processing jsj and quantum control of atomic and molecular processes Q . 
^ , Further, this analysis offers new insights into both decoherence and QCC, is relevant to semiclassical descriptions of 
decoherence dynamics [Q , and has motivated purely classical descriptions of dynamics-induced intrinsic decoherence, 
with preliminary computations Q that support the analytic results presented here. 

We begin by introducing classical analogs of some representation-independent and representation-dependent mea- 
^ . , sures of decoherence. A second-order perturbative treatment is then used to examine QCC in early-time decoherence 
' dynamics. Subsequently, we introduce a classical theory of strong decoherence that allows us to go beyond perturba- 
tion theory; the results are then compared with corresponding quantum theory. 

Consider a system with a time-independent Hamiltonian — P^/(2m) -I- V{Q), where {Q,P) are conjugate 
position and momentum variables, coupled to N independent harmonic bath modes described by the Hamiltonians 
Hj = p'j/{2mj) + rrijUjjq^ /2, where {pj,qj} (j = 1,2,---,A^) are bath-mode phase space variables. The system- 
bath coupling potential is assumed to be V^^ = ^jfiQ)lj that the total Hamiltonian is given by H = 
H" + T^^=i[^^j + Cjf{Q)lj\- The phase space distribution function evolved classically and the quantum Wigncr 
function for the entire phase space are represented by p'^[Q, P, {qj,pj}, t] and p^[Q, P, {qj,Pj},t], respectively. Their 
time evolution is given by dp'^/dt = {H, p^jp and dp^ /dt — {H, p^}m, where {■}p denotes classical Poisson bracket 
and {'Im denotes quantum Moyal bracket. Further, we define classical and quantum reduced distribution functions 
Pc{Q,P,t) = J p'^dT^ and pw{Q,P,t) = / p^dT^ , where dT^ = Ylj^i dqjdpj. Our interest is in the correspondence 
between the classical and quantum decoherence dynamics of initial states that can be either classical (i.e. positive 
Wigner density everywhere) or nonclassical [Q (e.g. displaying regions of negative p^[Q, P, {qj,pj}, t]). 

Consider two measures of decoherence in each of classical and quantum mechanics. One widely-used and 

representation- independent measure is the linear entropy Sq = 1 — Tr{p ) Q, where p is the reduced density operator 
of the system. An increase in Sq causes 1/(1 — 5*^) to increase, corresponding to an increasing number of incoherently 
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populated orthogonal quantum states. Since Sq = 1 — 27rh J p^{Q,P^t)dVs, where dVg = dQdP, this entropy has a 
natural classical analog (denoted Sc) obtained by replacing p^y with p^. That is, S'c = 1 — 2nh J pl{Q,P,t)dTg. A 
more detailed, but representation-dependent description of decoherence is the decay of off-diagonal density matrix 
elements such as {Qi\p{t)\Q2) ■ Significantly, we discover that the classical analog of these matrix elements can also 
be constructed. Specifically, noting that {Qi\p{t)\Q2) = J dPpw{Q, P, t) exp [iAQP/h], where Q = {Qi + Q2)/2 and 
AQ — Qi — Q2, we define the classical analog (denoted Pc(Qi, Q2, t)) of {Qi\p{t)\Q2) as the Fourier transformed clas- 
sical distribution function, i.e., Pc(QijQ2ji) = / dP pc{Q , P, t) exjp [iAQP/h]. This approach can be readily extended 
to the momentum representation. 

Perturbative treatments have proved to very useful in understanding decoherence dynamics [ p^ , 11 1. Here, to 



examine classical vs. quantum decoherence dynamics at short times, a regime of great interest in the control of 
decoherence, we consider a second-order expansion with respect to time variable t for both Sc and Sq, i.e., Sc{t) — 

Sc{0) + t/Tc,i + t^lTl_2 ^ , and Sq{t) ^ Sq{Q) + t/Tq^i + i^/r^ j H . Using the definitions of Poisson and Moyal 

brackets and assuming that the initial distribution function is decorrelated with initial bath statistics, we obtain 



4- = ^ / dgidg2|/5c(Qi,Q2,0)pAQ2 

'''c,2 " J 



df{Q) 



dQ 



(2) 



and 



4- = ^ / dQidQ2\{Qi\m\Q2)\^AQ^ 



A/(Q) 
AQ 



(3) 



where Cb = Ef=i C| coth(/3?iWj/2)/(2TOjWj), A/(Q) = f{Q + AQ/2)~f{Q~AQ/2), and /? is the Bohzmann factor. 



Note that the factor h appearing in the classical result [Eq. |^)] is just due to the definitions of Sc and pdQi, Q2, 0), 
and that the initial variances of the bath variables qj have been evaluated using quantum statistics to ensure the same 
initial quantum state for the ensuing classical and quantum dynamics. Note also that the decoherence time scale 
indicated in the easily-derived and simple quantum result of Eq. (^ is consistent with, but is more transparent than, 
a previous perturbation result (Eq. (5.6) in Ref. |12 ) obtained using a sophisticated influence functional approach. 

Equation (|^) shows that zero first-order decoherence rate i.e., 1/Tg^i=0, has a strict classical analog. More inter- 
estingly, Eqs. (^) and (||) show that, for the same fixed initial distribution function, the ratio of 1/T'g^2 to 1/t^2 
is ?i-independent. As seen from Eqs. @) and (^), (1/t^^2 ~ l/'^c.2) arises from the difference between the derivative 
df /dQ and the finite-difference function Af /AQ, weighted by AQ^ and the initial state. As a result: (1) For any given 
f{Q), as long as {Qi\p{0)\Q2) decays fast enough with \AQ\ such that Af / AQ « df /dQ, there would be excellent 
QCC in early-time decoherence dynamics. The smaller the Ti, the more rigorous is this requirement. (2) If f{Q) 
depends only linearly or quadratically upon the coupling coordinate Q, then (I/t^ 2 — l/''"c 2) — for any initial state. 
Significantly then, in all traditional decoherence models jl^ where f{Q) = Q is assumed, there exists perfect QCC in 
early decoherence dynamics, regardless of fi, and irrespective of the system potential V{Q) [|lj]. Indeed, in the case 
of f{Q) = Q Eq. (^ reduces to an important result, previously obtained quantum mechanically [l0| : 
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r^ = r^ = 2^>.:^3^coth( 



, I3hujj 



.7 = 1 ■' 



(4) 



where the initial state of the system is assumed to be pure, with the initial variance in Q given by 5^Q. (3) For 
nonlinear f{Q) where Af / AQ ^ df /dQ over the range of the initial state, QCC can be very poor. 

The second-order perturbative treatment is most reliable at short times and for weak decoherence. The results are 
particularly significant for studies of decoherence control where early-time dynamics of weak decoherence is important. 
In these circumstances it is useful to understand the extent to which (quantum) decoherence is equivalent to classical 
entropy production, i.e. to increasing Sc(t). In particular, if there exists good correspondence between classical and 
quantum decoherence dynamics, then the essence of decoherence control is equivalent to the suppression of classical 
entropy production, and various classical tools may be considered to achieve decoherence control. If not, then fully 
quantum tools are required. 

As an example, consider decoherence for an initial superposition state of two well-separated and strongly localized 
Gaussian wavepackets located at Qa = Qab ~ ^Qab/2 and Qh = Qab + ^Qab/2 with Q^f, = 0. For this initial state. 
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1/Tc,2 - {Cb/n)AQl, [df{Q,,)/dQ^,f, and l/r^^ ^ {a/n)AQl [A/(Q,,)/AQ,fcf . Then in a cubic decoherence 
model, for example, where f{Q) — , one would obtain ~ since d/((9afc)/'^'9ab = 0- However, here 

^/Tq2 l/'''c 2' there is appreciable decoherence without classical entropy production. By contrast, in another 
nonlinear decoherence model where f{Q) — sin{2'KQ/ AQab + 7r/4), l/r^2 ~ since f{Qa) — f{Qb)- Here, however, 
^I'^c2 l/'''i3 2j the system is decoherence- free hut with substantial classical entropy production. Since we find 
that the ratio of r^j to t^2 early-time decoherence dynamics is independent of % for fixed initial state, these two 
examples lead to a rather counter-intuitive result: given a macroscopic object which is initially in a superposition 
state of two distinguishable states and is nonlinearly coupled with an environment, classical dynamics could totally 
fail to predict its initial entropy production or its decoherence rate. Indeed, Eqs. ^ and (^) suggest that, as 
long as df{Q)/dQ ^ and |/(Q)| is bounded, then 1/t^2 saturates with increasing AQab, whereas 1/t^2 does not. 
Thus, one can conclude that decoherence dynamics must be quantum and that the system-environment coupling 
must be nonlinear if the saturation behavior of early-time decoherence rates is observed experimentally [p^ . Further, 
it is clear that in the limit of large AQab, classical decoherence dynamics in the general case of nonlinear system- 
environment coupling predicts much faster decoherence than does quantum decoherence dynamics. This leads to the 
rather surprising inference that initial superposition states of well-separated wavepackets would be more susceptible 
to nonlinear system-environment coupling if they are propagated by classical dynamics than by quantum mechanics. 

To go beyond the perturbation results we now consider a strong decoherence model in which decoherence is assumed 
to be much faster than the system dynamics, so that can be set to zero |Q. We consider both the "off-diagonal 
elements" pc{Qi,Q2,t) as well as the entropy Sc(t) and compare them to the quantum results. 

In this case the classical Liouville dynamics gives 

9F,[Q, AQ, t] dHl dF,\Q, AQ, t] 

dt f-; dqk dpk 

^ dHl 9F,[Q,AQ,fa-,p,},t] ^ dF^^QMi^P^lA 
^ opk dqk ^ dpk 

-^AQj2Ck^^qkF4Q,AQ,{q„p,},tl (5) 

where Fc{Q, AQ, {qj,pj},t) = J dP exp[iAQP/h]pc[Q, P, {qj,Pj},t]. Since Q — due to Hs = 0, and AQ is a 
time-independent parameter introduced in the Fourier transformation, Eq. (||) leads to 

dF,[Q, AQ,{qjit),pj{t)},t] ^ dF,[Q, AQ,{q,{t),p,{t)},t] 

dt dt 
, ^ dF,[Q,AQ,{qk{t),pk{t)},t] . , ^ dF,[Q,AQ,{qk{t),pk{t)},t] . 

k^l ^^f"^*' k^l ^P''^^) 

= -jAg^Cfegfc(t)^F,[Q,AQ,{g,(<),p,(i)},t], (6) 

IT- k=l 

where {qj{t),pj{t)} satisfy qj{t) ~ dHj /dpj{t) and Pj{t) — —dH^ /dqj{t) — Cjf{Q), of which the solution is 

qj (t) = '^^■^'^^J [cos{ujjt) - 1] + qj (0) cos(cjji) + sm{ujjt), (7) 

and Pj{t) = mjqj{t). Analytically integrating Eq. (^, and using dV^ (t) = dV^ (0) and pc{Qi,Q2,t) = 
/ dT^{t)Fc^, AQ, {qjit),pj{t)}, t], we have 

Pc{Qi,Q2,t)^ J dT^{0)F,[Q,AQ,{qjiO),p,iO)},0] 

xeM-^T f dtAQY^Ck'^qkit)], (8) 
Substituting Eq. (0) into Eq. (|^), using the initial quantum state of the bath that is initially uncorrelated with the 
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system, and assuming that the equihbrium state of the bath is maintained, we obtain 



Pc{Qi,Q2,t) 
/5c(Qi, Q2,0) 



exp 



where ^t) = {AQ)fiQ)[df{Q)/dQ]Bi{t)/h, with Bi{t) = E1i ^^H* - sin^OMl/C 



C| coth(/3?iWj72)[l — cos(ajjt)]/(2mj?iw|). Interestingly, the classical result [Eq. (^] displays two dynami- 
cal aspects of pc{Qi, Q2,t), i.e., coherent dynamics of its phase (l)c(t), and incoherent decay due to bath statistics. 
The classical linear entropy Sc{t) can then be obtained from Eq. as 



(9) 



and B2{t) = 



Sc{t) = 1-/ dQidQ2\pciQi,Q2,0)\^e^p 



With similar manipulations for quantum strong decoherence dynamics, we obtain the quantum result 



{Qi\mQ2) 

(Qi|p(0)|Q2) 



exp 



where = AQf{Q)[Af{Q)/AQ]Bi{t)]/n, and 



Sg{t) = 1 - / dQidQ2|(Qi|p(0)|Q2)pexp 



^2(i) 



(10) 



(11) 



(12) 



These results extend those in Ref. to nonlinear f{Q) using a simple approach and demonstrate a direct classical 
analog to quantum strong decoherence dynamics. 

Since dB2{t)/dt{t = 0) = and d'^ B2{t) / dt^ [t = 0) = Cb/h, one finds that in the short time hmit, Eqs. ([lO|) and 
( p^ ) reduce to previous perturbation results of l/rc,i, l/r^2: l/'^q.ij and 1/t^2- Furthermore, the classical results 
[Eqs. (|) and (|o|)] are again much similar to the quantum results [Eqs. ( pi] ) and dl^)], with the only difference being 
that A//AQ in the quantum expression is replaced by df /dQ in the classical result. 

This result makes clear that our previous QCC results based upon second-order perturbation theory 
are generalizable to all orders of time in the strong decoherence case. In particular, defining 7c (i) = 

dln|pc((3i,Q2,t)|M and 7,(t) = din |(Qi|^(t)|Q2)lM, we have 7^(0 = -{AQf [df(^)/dQ\' {dB2{t)/dt), and 
lq{t) = —{AQY [A/((5)/A(5] {dB2{t)/dt). Then, in the case of hnear and/or quadratic coupling, e.g., f{Q) = 
aQ + bQ^, one has 7c (i) = 7q{t) and Sc{t) = Sq{t), showing that there is perfect QCC in decoherence dynamics for 
all times. 

By contrast, in the case of nonlinear coupling, "fdt) in general does not saturate with increasing AQ whereas 
jq{t) does saturate for bounded \f{Q)\- As such, in the limit of large AQ, one has |7c(i)| >> hq{t)\ and thus 
[1 — Sc{t)] << [1 — Sq{t)] as t increases, with \jc{t)/'jq{t)\ independent of h. This observation is of conceptual 
importance: it says that decoherence can dramatically improve QCC, but as far as some detailed characteristics 
of decoherence dynamics are concerned, decoherence itself does not necessarily suffice to ensure that the dynamics 
of quantum entropy production equals that of classical entropy production. That is, even in the presence of strong 
decoherence, subtle quantum classical differences may persist in some measures (e.g., 1/[1 — S'g(t)] vs. — Sc{t)]) for 
all finite times. Note, however, the entropy measures such as 1/[1 — Sq{t)] are not a quantum mechanical observables 
and hence do not allow one to directly measure the subtle difference between classical and quantum decoherence 
dynamics at later times. 

Thus, from both the perturbation and strong decoherence results, we obtain that QCC depends critically upon the 
initial quantum state and the nature of the system-environment coupling. This result should have an impact on our 
current understanding of decoherence even when the role of the dynamics of the system is important. For example, it 
is worthwhile reexamining the relationship between classical Lyapunov exponents and decoherence rates in classically 
chaotic systems, since previous studies [O] only dealt with the case of linear system-environment coupling. 

In conclusion, we have examined, using analogous measures, the decoherence dynamics of an initial quantum state 
coupled to a bath that is subjected to either classical or quantum dynamics. Within the framework of a second- 
order perturbative treatment and a strong decoherence theory, we have exposed the system-independent conditions 
under which the quantum decoherence dynamics is either well, or poorly, approximated by classical dynamics. Further 
studies are ongoing to assess QCC in cases beyond the short time and strong decoherence approximations. Preliminary 
computational results |^ support the conclusions drawn herein. 

This work was supported by the U.S. Office of Naval Research and the Natural Sciences and Engineering Research 
Council of Canada. 
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